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We discuss a novel method to calculate /s on the lattice, introduced in Q, based on the study of the dependence 
of finite size effects upon the heavy quark mass of flavoured mesons and on a non-perturbative recursive finite size 
technique. This method avoids the systematic errors related to extrapolations from the static limit or to the tuning 
of the coefficients of effective Lagrangian and the results admit an extrapolation to the continuum limit. We show 
the results of a first estimate at finite lattice spacing, but close to the continuum limit, giving /a = 170(11)(5)(22) 
MeV. We also obtain /s^ = 192(9)(5)(24)MeV. The first error is statistical, the second is our estimate of the 
systematic error from the method and the third the systematic error from the specific approximations adopted in 
this first exploratory calculation. The method can be generalized to two-scale problems in lattice QCD. 



1. Introduction 

Lattice QCD evaluations of quantities charac- 
terised by two scales with a large hierarchy re- 
quire in general a very high lattice resolution and 
a sizeable total physical volume to correctly ac- 
count the dynamics of the small distance scale 
and to dispose of the finite size effects related to 
the large distance scale. A good example is pro- 
vided by the pseudoscalar B meson decay con- 
stant [^, where the small distance scale is repre- 
sented by the inverse of the bottom quark mass 
and the large distance scale by the radius of the 
B meson, related in turn to the inverse of the 
light quark mass. A straight evaluation of the de- 
cay constant would require lattices with N — 80'* 
points or more, exceeding the present generation 
computers capabilities, and, in the case of un- 
quenched simulations, the ones of the next gen- 
eration. One resorts to approximate calculations 
based on extrapolations from the static limit or on 
non-relativistic formulations of standard QCD. 
All the available methods introduce systematic 
errors related to extrapolation fits and/or to the 
use of effective Lagrangians. 

We discuss a novel approach based on the study 
of the dependence upon the heavy quark mass of 
finite size effects for the pseudoscalar decay con- 
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stant of heavy flavoured mesons (see [Q for de- 
tails). The basic assumption is that the finite 
size effects are mainly related to the light quark 
mass and rather insensitive to the one of a suffi- 
ciently heavy quark. We discuss the general fea- 
tures of the method assuming the continuum limit 
has been taken. The relevant quantity is the ratio 
cr = f b{'2,L) / f b{L) of the pseudoscalar constants 
at different volumes, where fsiL) is the value of 
the decay constant on a volume with linear size L. 
The dimensionless a depends on general grounds 
upon three dimensionless variables: miL, rrihL 
and AqcdL. For a sufhciently large heavy quark 
mass m/j, the dependence is basically dominated 
by the light quark and the expansion for large rrih 
takes the form 

t T A ^^ , C (meL, AqcdL) 
a = a [miL,AQCDL) -\ . (1) 

A simple phenomenological ansatz for a can be 
made based on the concept of a reduced mass con- 
structed out of the heavy and light quark masses 

a = a (niredL, AqcdL) , (2) 

where rrired — (a*iM2)/(mi + fJ-2)- The quantity 
fj,i is a function of the quark mass, but not only: 
indeed, for very light masses, finite size effects 
are regulated by the physical meson size, which 
is expected to remain finite when the light quark 
mass approches zero. 
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A crucial question is the threshold value of the 
quark mass on a given volume where the large 
mh expansion becomes reliable. As we will show, 
this value falls in a mass range of the order of 
a couple of GeV in the renormalization invariant 
mass scheme, where the calculation on a single 
lattice is affordable. Under these circumstances, 
the strategy to obtain /s is the following. One 
first performs a calculation on a lattice where the 
resolution is suitable for b quark propagation, but 
the total volume is unavoidably a small one. This 
sets /b on a finite volume. In order to connect to 
the large volume results, one needs the step scal- 
ing function a for values of heavy quark masses 
generally lower than those of the simulation where 
the finite size value of fs was obtained. The 
possibility of extrapolating a to heavier masses 
depends upon the validity of the asymptotic ex- 
pansion: in a favourable case, as is the real one, 
one can evaluate the finite size effects in a reli- 
able way, connecting, by a repeated iteration of 
the procedure, small volume values of /b to the 
ones on large volumes, 

rJ''^'^fB{LoMLoM2Lo)... , (3) 

and the recursion stops on a volume where 
CT ~ 1 within a required precision. The con- 
tinuum limit is obtained by extrapolating to 
zero lattice spacing the step scaling function ob- 
tained at fixed physical quantities. We use the 
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Figure 1. Step scaling function Eo.4-0.8 for the evo- 
lution of fhe from 0.4 fm to 0.8 fm at /3 = 6.737. 

Schrodinger Functional framewrok to compute 
the pseudoscalar decay constants and we fix the 



lattice scale in terms of ro Meson masses, on 

different volumes, are tuned by fixing the physi- 
cal value of the Renormalization Group Invariant 
quak masses 

2. Results and Discussion 

The results are obtained at finite lattice spac- 
ing. The size of the smallest volume follows from 
the decision of making our estimate for the finite 
size Jb on a 48 X lattice with a cutoff of about 

^ ~ 12 GeV. The value of the bare couphng 
for this lattice spacing has been obtained from a 
fit in ref. |Q. The procedure fixes /3(ao) = 7.3 and 
the physical volume Lq = 0.4 fm. On this lattice, 
we simulate heavy quark masses up to 0.3 in lat- 
tice units, corresponding to bare physical masses 
slightly above 4 GeV. Indeed, as a general caution 
against large lattice artifacts, at all P values we 
take the maximum heavy quark mass in lattice 
units of the order of 0.3. The first E (we dis- 
tinguish between the continuum step function a 
and the one at finite lattice spacing E) goes from 
the volume of 0.4 fm to the one of 0.8 fm. In 
terms of lattice points, we go from 12 to 24, and 
we have to match the starting volume of 0.4 fm 
with a resolution which is half of the one used for 
a correct estimate of the bottom quark propaga- 
tion. According to our caveat, it follows that the 
maximum bare quark mass that we can achieve 
is correspondingly halved, i.e. of about a cou- 
ple of GeV at a bare couphng /3 = 6.737. We 
make a further iteration with a second E going 
from 0.8 fm to 1.6 fm, where our investigation 
of heavy quark masses stops at the order of the 
charm quark mass. The corresponding bare cou- 
pling is (3 = 6.211. The finite volume effects for 
this second evolution step are small enough to 
make the neglection of the residual volume effects 
a safe assumption, that however can be tested ex- 
plicitly. 

The plots in Figs. 1 and 2 show the dependence 
of E upon the heavy RGI quark mass M^qj for 
the two volume jumps and provide evidence for 
a plateau of insensitivity to heavy quark masses: 
the data have been obtained from a linear extrap- 
olation in M^Qj to the down and strange RGI 
quark masses reported in 0. These figures, with 
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the figures of ref.||l|], support the procedure pro- 
posed. The statistical errors are computed by a 
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Figure 2. Step scaling function Eo.g-i.e for the evo- 
lution of fhi from 0.8 fm to 1.6 fm a,t f3 — 6.211. 

jacknife method (see ref. |^] for details). The fi- 
nite size value of Jb is obtained by a calculation 
on the highest resolution lattice and the RGI bot- 
tom quark mass is obtained from the equation 



M = ZMi9Q)mwii9o) 



(4) 



In order to obtain the renormalisation constant 
Znigo) at /3 = 7.3 and P — 6.737, we have used a 
safe interpolation of the pseudoscalar renormal- 
isation constant Zp{go,^) at a value of fi three 
times the reference value used in eq. (6.8) of 
ref. M. The value for /b that we obtain is 



/i3(0.4fm) = 483(4)MeV 



(5) 



By using the values of S for the b quark at con- 
stant RGI mass 



^0.4- 



0.401(4), 



0.88(4) (6) 



we obtain our estimate of /b on the large volume: 



phys 

B — 



170(ll)MeV 



(7) 



where the error quoted in the previous equation 
is statistical only. In this the written version of 
the talk we include our estimate of the sistem- 
atic errors that can be partly ascribed to specific 
approximations used in the present computation 
that can be eventually removed, and partly to 
the uncertainty in the extrapolation in the heavy 
quark mass of finite size effects, inherent to the 



method proposed. The overall error on the num- 
ber Jb coming from the removable systematic un- 
certainties is of about 13% and of at most 2 — 3% 
from the ones deriving from the unavoidable ex- 
trapolation in the heavy quark mass, leading to 
a global uncertainty of about 25 MeV of which 
about 20 are removable while 5 stay with the 
method: 



fphys 
J B 



170(ll)(5)(22)MeV . (8) 
In the same way we obtain 
jj^ys ^ i92(9)(5)(24)MeV . (9) 
We quote also the results of the ratios 



1.13(2)(1) 



rphys 
J D, 

fphys 
J D 



1.10(1)(1) (10) 



rphys 
J B, 
rphys 
J B 

that were asked at the end of the talk. 

The method proposed can be generalized to 
problems characterised by two very different mass 
scales, if the decoupling of the large mass scale 
from the low scales of non-perturbative QCD dy- 
namics holds true. This appears to be the case in 
the example discussed and is somehow supported 
by the wide success of the predictions of perturba- 
tive QCD calculations for hard processes that are 
insensitive to the dressing mechanism of quarks 
and gluons into standard hadronic final states. 
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